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Abstract 

Standard analyses of low-energy NN and nuclear parity- violating observables 
have been based on a n— , p—, and a;— exchange model capable of describing 
all five independent s — p partial waves. Here a parallel analysis is performed 
for the one-body exchange-current, and nuclear polarization contributions to 
the anapole moments of 133 Cs and 205 T1. The resulting constraints are not 
consistent, though there remains some degree of uncertainty in the nuclear 
structure analysis of the atomic moments. 



1 



Nuclei and nucleon-nucleon scattering are the only experimentally tractable systems in which 
to study the flavor- conserving hadronic weak interaction, where neutral current effects arise. 
This interaction can be isolated, despite the presence of much larger strong and electromag- 
netic effects, because of the accompanying parity violation. The long-term goal of the field 
is to learn how standard-model quark-boson couplings give rise to long-range weak forces 
between nucleons |l]-|3[]. 

Several precise and interpretable measurements of parity nonconservation (PNC) in nu- 
clear systems have been made. These include the longitudinal analyzing power A z for p + p 
at 13.6 and 45 MeV, A z for p+a at 46 MeV, the circular polarization P 7 of the 7-ray emitted 
from the 1081 keV state in 18 F, and Ay for the decay of the 110 keV state in polarized 19 F. 
An analysis of these results, which have been in hand for some time, suggests that the 
isoscalar PNC NN interaction is comparable to or somewhat stronger than the "best value" 
suggested theoretically, while the isovector PNC NN interaction is significantly weaker, an 
isospin anomaly superficially reminiscent of the AI = 1/2 rule in strangeness-changing de- 
cays. There has been great interest in obtaining additional experimental constraints that 
would test this tentative conclusion, as the weakness of the isovector weak NN interac- 
tion was unexpected. It had been anticipated that the neutral current would enhance this 
component (which is dominated by 7r exchange). 

After a considerable wait, several new PNC measurements have become available or are 
expected soon. Recently the Colorado group f| measured, for the first time, a nuclear 
anapole moment - the PNC axial coupling of a photon to the nucleus in its ground state 
- by determining the hyperfine dependence of atomic parity violation. A significant limit 
on the anapole moment of another nucleus, 205 T1, has also been obtained p[. Preliminary 
results || for the p + p A z (221 MeV) are now available: this result, which is sensitive only to 
PNC p exchange, can be combined with lower energy measurements to separately determine 
the p and uj PNC couplings. In the next few years new results are expected on the PNC 
spin rotation of polarized slow neutrons in liquid helium M and oni 7 in n + p — > d + 7 || . 

The primary obstacle to an analysis in which the new PNC constraints are combined 
with older results is the difficulty of treating anapole moments with a comparable degree 
of sophistication. The theoretical framework for NN and nuclear observables is a it—, p— , 
and uj— exchange model involving six weak meson-nucleon couplings f n , h° p , h 1 , h 2 p , h^, and 
/i*. (The superscripts denote the isospin.) We use the standard Desplanques, Donoghue, 
and Holstein (DDH) definitions of these couplings [|IJ (so the sign of f n differs from that 
of Ref. 0). If one takes the view that the PNC observables involve a momentum scale 
comparable to or below the pion mass, then this framework is quite general, describing 
the five independent s — p amplitudes and the separate long-range 7r contribution to those 
amplitudes. 

The older PNC results involve systems that are either amenable to exact potential-model 
calculations, or can be "calibrated" by a relation between PNC mixing and axial charge (3 
decay @. Thus it is generally believed that little nuclear structure uncertainty remains. 
Calculations account for the full two-body PNC potential and the effects of short-range 
correlations on the potential. In contrast, most anapole moment investigations have been 
evaluated in the extreme single-particle (s.p.) limit employing effective one-body potentials. 
The only calculation employing a modern strong effective interaction (in a large-basis shell 
model (SM) calculation) in combination with the PNC two-body potential was limited to 
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the effects of f n 0. Here that calculation (for 133 Cs) is extended to the full potential and 
then repeated for 205 T1. We then examine the consistency of these and other constraints on 
the weak potential. 

Zel'dovich [|10J introduced the anapole moment, the El coupling of a virtual photon to 
an elementary particle induced by PNC. The measurability of anapole moments was first ar- 
gued by Flambaum, Khriplovich, and Sushkov JTI] , who studied the nuclear-spin-dependent 
vector(electron)-axial(nucleus) interaction induced by the nuclear anapole moment in atomic 
PNC experiments. A definitive extraction of the anapole contribution to atomic PNC eluded 
experimentalists until the 7a result of Wood et al. on 133 Cs 0. 

Electrons in an atom experience a weak contact interaction with the nucleus of the form 

G -> 

H w = -JLna- Ip(r), (1) 

— # 

where / and p(r) are the nuclear spin and density. (Note that k differs from the definition 
of HITi[l2|.) From the hyperfine dependence of the atomic PNC signals in 133 Cs (as extracted 
by Flambaum and Murray JT2|) and 205 T1 one finds 

K ( 133 Cs) = 0.112 ±0.016 

k( 205 T1) = 0.293 ± 0.400. (2) 
One contribution to k originates from Zq exchange with axial coupling to the nucleus 

2 v WJ (/||/||/> 

where qa = 1-26 is the axial vector coupling, sin 2 9w = 0.223, and || denotes a matrix ele- 
ment reduced in angular momentum. The reduced matrix element of I is yl{l + l)(2i + 1). 
The Gamow- Teller matrix elements, taken from the SM studies described below, are -2.305 
( 133 Cs) and 2.282 ( 205 T1), not too different from the corresponding s.p. values of -2.494 (un- 
paired l<77/2 proton) and 2.449 (3si/ 2 proton). Thus the predicted k Zo are 0.0140 and -0.127, 
respectively. Note that one-loop standard model electroweak radiative corrections modify 



this result, reducing the isovector contribution and inducing a small isoscalar amplitude |L3| 



A second contribution to k is generated by the combined effects of the usual coherent 
Zq coupling to the nucleus (vector coupling, proportional to the nuclear weak charge Qw) 
and the magnetic hyperfine interaction ||14|| . From the measured nuclear weak charge and 



magnetic moment Bouchiat and Piketty [15] find 



^( 133 Cs) = 0.0078 

^( 205 T1) = 0.044. (4) 

Thus the experimental values for the anapole contributions to k are obtained by subtracting 
the results of Eqs. (3) and (4) from Eq. (2), yielding 

KanapoleC^Cs) = 0.090 ± 0.016 

Tl) = 0.376 ± 0.400. (5) 



K 



anapolei 



(205 
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These values can then be related to the corresponding nuclear anapole moments by 



47ra v / 2(I||i 1 ||/)/e 



■±iuxyj, \-t ||^i||J//c , > 

^anapole ^/f2 n ir\\r\\r\ ' V / 



where v4ia is the anapole operator. As discussed in Ref. ||, the anapole moment is the 
sum of three contributions, the one-body operator associated with the anapole moments 
of individual nucleons, the two-body exchange currents, and the polarization contribution 
associated with parity admixtures in the nuclear ground state. The goal is to treat each of 
these in a manner consistent with the standard DDH meson-exchange model of PNC. 

The anapole operator can be written, via the extended Siegert's theorem, in a form where 
all components of the current that are constrained by current conservation are explicitly 
removed. This yields 



Jtx(r) + v 7 ^ (Y 2 (n r ) ® jf» (#0) J • (7) 



(This form is equivalent to others in common use when current is conserved, but this is 
generally only possible in the simplest nuclear models, such as the s.p. limit.) This El 
operator has nonzero ground state matrix elements for currents induced by PNC, as well as 
for the ordinary electromagnetic current through PNC admixtures in the ground state. 

a) Nucleon anapole moment. The one-body PNC electromagnetic current is obtained from, 
for example, loop diagrams involving a '-fim vertex, where one pion is absorbed by the 
nucleon with a PNC coupling and the other with a strong coupling, yielding 

A nucleonic = £ ^ + ^(Q)^)] ^(j). ( 8 ) 
i=l 

It is thus clear that the contribution of spin-paired core nucleons cancel, leaving only the 
anapole moment of the valence nucleons. In our earlier work |J only the pion contribution 
to a s (0) and a v (0) was included, yielding a result proportional to ef-^g^NN, where g n NN is 
the strong coupling. The isoscalar coupling a s (0) proved to be about four times larger than 
a v (0). This work was recently extended to included the full set of one-loop contributions 
involving the DDH vector meson PNC couplings, using the framework of heavy baryon chiral 
perturbation theory and retaining contributions through 0(1/A^), where A x = A7tF n ~ 1 
GeV is the scale of chiral symmetry breaking |Tj| . The contributions due to f n are consistent 
with the earlier work - the new a s (0) is about 1.3 times larger, while a v (0) is zero in this 
order - while the addition of the heavy mesons greatly enhances a v (0). An evaluation with 
DDH best value couplings yields a v (0) ~ 7a s (0). Folding the resulting expressions with our 
SM matrix elements ((/|| Y*Li — -2.372 and 2.532 for Cs and Tl, respectively) yields 

the results in Table I. The inclusion of heavy meson contributions substantially enhances 
the one-body anapole terms and alters the isospin character, generating opposite signs for 
the proton and neutron anapole moments. 

b) Exchange currents. Insertion of the NN pair and transition currents, where the meson 
exchange involves a PNC coupling on one nucleon and a strong coupling to the second, 
into Eq. (7) produces a two-body PNC anapole current. The only previous estimate || 
of contributions of this type was restricted to pions. The extension to include the p and uj 
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TABLE I. Decomposition of the SM estimates of the anapole matrix element (I||Ai||I)/e into 
its weak coupling contributions. 



Nucleus 


Source 


fn 






hi 


hi 


hi 


133 Cs 


nucleonic 


0.59 


0.87 


0.90 


0.36 


0.28 


0.29 




ex. cur. 


8.58 


0.02 


0.11 


0.06 


-0.57 


-0.57 




polariz. 


51.57 


-16.67 


-4.88 


-0.06 


-9.79 


-4.59 




total 


60.74 


-15.78 


-3.87 


0.36 


-10.09 


-4.87 


205 rpj 


nucleonic 


-0.63 


-0.86 


-0.96 


-0.35 


-0.29 


-0.29 




ex. cur. 


-3.54 


-0.01 


-0.06 


-0.03 


0.28 


0.28 




polariz. 


-13.86 


4.63 


1.34 


0.08 


2.77 


1.27 




total 


-18.03 


3.76 


0.33 


-0.30 


2.76 


1.26 



PNC couplings is a formidable task requiring evaluation of the p and u pair currents and 
the pfrf and ptfy currents. The first step, to identify which of these currents are important, 
was a Fermi gas reduction of the operators to one-body form, which requires a standard 
summation over occupied core states. The results indicated that the ppj, pnj, and the 
component of the u pair current where the photon and PNC uj couplings are on different 
nucleon legs are negligible, well below 1% of the dominant tt currents. We then evaluated 
the remaining (and more important) terms using the two-body density matrices from our 
large-basis SM calculations. Finally, to cross check our procedure for the neglected currents, 
we compared the ratios of two-body results to the ratios of the Fermi-gas reductions of the 
same currents. The correspondence was typically at the few percent level, which supports 
the notion that the one-body approximation could be used to assess the relative importance 
of all currents. (The approximation is less useful as an estimator of absolute strengths [0.) 
The exchange current totals are given in Table I. It is clear that the 7r contribution continues 



to dominate. This work in described in considerable detail in Ref. |16 |. 

c ) Nuclear polarization contribution. The nuclear polarization contribution to the anapole 

moment is given by 

(I\\Ar\\n){n\H PNC \I) + ^ 

where A^ 71 is obtained from the ordinary electromagnetic current operator, \I) is a ground 
state of good parity, H PNC is the PNC NN interaction, and the sum extends over a complete 
set of nuclear states n of angular momentum I and opposite parity. Several nontrivial nuclear 
structure issues must be addressed in evaluating this term. 

One of these affects all of our calculations, the quality of the ground state wave function. 
The canonical SM space for 133 Cs is that between the magic shells 50 and 82, L77/2 — 20I5/2 — 
l^n/2 — 3si/2 — 2g?3/2. Calculations were performed with protons restricted to the first two 
of these shells and neutron holes to the last three, producing an m-scheme basis of about 
200,000. Two interactions were employed, the Baldrige-Vary potential employed in Ref. || 
and a recent potential developed by the Strasbourg group ]T7), both of which are based 
on the addition of multipole terms to g-matrix interactions and are designed for the 132 Sn 



region. As the results are very similar [16 , we quote only the former here. 2U5 T1 is described 
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as a proton hole in the orbits immediate below the Z=82 closed shell (3si/ 2 — 2d 3 / 2 — 2d 5 / 2 ) 
coupled to two neutron holes in valence neutron space between magic numbers 126 and 82 
(3pi/2 — 2/5/2 — 3p 3 / 2 — U13/2 — 2/7/2 — 1^9/2)- A Serber- Yukawa force was diagonalized in 
this space. 

The summation over a complete set of intermediate states in such spaces is impractical 
either directly or by the summation-of-moments method discussed in Ref. ||. However, 
because no nonzero El transitions exist among the valence orbitals, an alternative of com- 
pleting the sum by closure, after replacing 1/E n by an average value (1/E), is quite attrac- 
tive: the resulting product of the one-body operator Af 71 and two-body H PNC contracts to 
a two-body operator, so that only the two-body ground state density matrix is needed. 

The closure approximation can be considered an identity, clearly, if one knows the correct 
(1/E), e.g., if it were governed by nuclear systematics and thus could be predicted reliable. 
This is known to be the case for the various sum rules of the El operator, an operator with 
very similar properties to A e ^ and to the one-body analog of H PNC , a ■ p. To assess the 
situation for the anapole polarizability, we completed a series of exact calculations in lp— 
and light 2sl<i— shell nuclei ( 7 Li, 11 B, 17 > 19 > 21 F, 21,23 Na), determining the ground states from 
full Ohu diagonalizations. After performing the summations (by Lanczos moments methods 
||) over the lhu spaces, the dimensions of which range up to ~ 0.5M, the results in Table II 
were obtained. The systematics of the relationship between the closure results and the 1/E- 
weighted sums are both interesting and encouraging. Distinct excitation energies are needed 
for the three isospins contributing to H PNC . Measured as a fraction of the 1/E- weighted 
giant dipole average excitation energy, which is (1/E)^ 1 ~ (22-26) MeV for these nuclei, 
the appropriate effective energies for the closure approximation are 0.604 ± 0.056 (h° p , h^), 
0.899 ± 0.090 (/„.), and 1.28 ± 0.14 (h 2 p ). The larger (1/E) for h° and hi enhances these 
contributions to the anapole polarizability. The small variation in (1/E), once the isospin 
dependence is recognized, supports the notion that we can connected the closure result to 
the true polarization sum. From the known El distribution |T^] in 133 Cs we then determine 
T=0,l,2 closure energies of 9.5, 14.1, and 20.2 MeV, respectively. That is, we fix these as 0.6, 
0.9, and 1.28 of the El closure energy evaluated from the experimental dipole distribution. 
The corresponding 205 T1 values are 8.7, 12.9, and 18.5 MeV. The resulting polarization 
contributions are given in Table I. 

A summary of PNC constraints is presented in Table III and Fig. 1. Although the PNC 
parameter space is six-dimensional, two coupling constant combinations, — 0.12/i^, — 0.18/^^ 
and h° n + 0.7h° n dominate the observables, as Table III illustrates. We include the results for 

P tut ' 

A p z p at 13.6 and 45 MeV and the soon-to-be-published results at 221 MeV, A p z a at 46 MeV, 
P 7 ( 18 F), A),( 19 F), and the Cs and Tl anapole results. We do not include P 7 ( 21 Ne) because 
of the arguments of Ref. that the underlying nuclear structure is uncertain. The la error 
bands of Fig. 1 are generated from the experimental uncertainties, broadened somewhat 
by allowing uncorrelated variations in the parameters in the last four columns of Table III 
over the DDH broad "reasonable ranges." The resulting pattern is disconcerting. Before the 
anapole results are included, the indicated solution is a small /„- and an isoscalar coupling 
somewhat larger, but consistent with, the DDH best value, — (h° p + 0.7hl) DDHb ' v ' ~ 12.7. 
The anapole results agree poorly with the indicated solution, as well as with each other. 
Although the Tl measurement is consistent with zero, it favors a positive anapole moment, 
while the theory prediction is decidedly negative, given existing PNC constraints. The Cs 



6 



TABLE II. Functional dependence of SM results for the unweighted and 1/E- weighted anapole 
polarization sums and extracted T=0,1 average excitation energies, measured as a fraction of the 
{1/E) derived from the corresponding El polarizability sums. These are odd-proton nuclei. Note 
the closure result faithfully reproduces the correct /i° — /i° combination, but that distinct T=0,1 
excitation energies must be used with the closure result to reproduce the correct 1/E- weighted 
isoscalar/isovector ratio. In the case of 19 F, the lowest, nearly degenerate 1/2 - state was removed 
from all sums. 



Nucleus Unweighted 1/E- weighted (V-E)t=o ( 1 / e )t= 

^Li f n - 0.25(/i° + 0.63/i° ) U ~ 0.34(/i° + 0.58/i° ) 059 0.80 

n B f w - 0.42(/»° + 0.52/i° ) f w ~ 0.53(/»° + 0.52/i°) 0.70 0.89 

17 F /„.- 0.39(/i° + 0.46/i°) / 7r -0.60(/i° + 0.48/i°) 0.66 1.02 

19 F / 7r -0.21(/t° + 0.59/i°) / w -0.33(/i° + 0.56/i°) 0.58 0.90 

21 F / w - 0.25(/»° + 0.54/i° ) / w -0.41(/i° + 0.55/i°) 0.60 0.97 

21 Na / w - 0.40(/i° + 0.49/i° ) f n ~ 0.57(/i° + 0.51/i° ) 0.54 0.77 

23 Na / w - 0.40(/i° + 0.52/i° ) f n - 0.67(/i° + 0.53/i° ) 0.57 0.95 



TABLE III. PNC observables and corresponding theoretical predictions, decomposed into the 
designated weak-coupling combinations. 



Observable 


Exp.(xl0 7 ) 


U-0.12hl-0.18hi 


h° p + 0.7/£ 


h p 




h° 


hi 


APP (13.6) 


-0.93 ± 0.21 




0.043 


0.043 


0.017 


0.009 


0.039 


A»(45) 


-1.57 ± 0.23 




0.079 


0.079 


0.032 


0.018 


0.073 


Af> (221) 


prelim. 




-0.030 


-0.030 


-0.012 


0.021 




AP Q (46) 


-3.34 ± 0.93 


-0.340 


0.140 


0.006 




-0.039 


-0.002 


P 7 ("F) 


1200 ± 3860 


4385 




34 






-44 


^4 ( 19 F) 


-740 ± 190 


-94.2 


34.1 


-1.1 




-4.5 


-0.1 


(\\Ai\\)/e, Cs 


800 ± 140 


60.7 


-15.8 


3.4 


0.4 


1.0 


6.1 


(PalD/e, Tl 


370 ± 390 


-18.0 


3.8 


-1.8 


-0.3 


0.1 


-2.0 



result tests a combination of PNC couplings quite similar to those measured in A 7 ( 19 F) and 
in A? 01 , but favors larger values. 

The nature of the anapole discrepancy - theory constrained by other PNC measurements 
requiring a smaller Cs moment - is surprising. The first criticism of the theory would be that 
the SM calculations are still too limited, not generating the proper quenching of operators 
like (XT3 that are known to be sensitive to core polarization effects. There is evidence, in 
the case of Tl where the odd proton is identified with the I = 3si/ 2 orbital, that this is 
the case: the SM predicts a spin-dominated magnetic moment of 2.59 fiN, improved from 



the s.p. value 2.79 fiN but well above the experimental value 1.64 /i^r. Indeed, in Ref. (15 
similar arguments were used to invoke quenching factors for the s.p. anapole moment results. 
We have explored these issues using the light nuclei of Table II as test cases, and while the 
approach of Ref. |L5| proves somewhat naive, the guess about the quenching effects of core 



polarization is qualitatively correct |16] . Clearly if anapole matrix elements are further 
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- 0.12 h p - 0.18 h w 

FIG. 1. Constraints on the PNC meson couplings (xlO 7 ) that follow from the results in Table 
III. The error bands are one standard deviation. 



quenched, the Cs problem in Fig. 1 becomes worse. Anticipated new results on np PNC - 
from the n+ 4 He and n + p — > d + 7 experiments - will be of great help. 

Our numerical results for Cs are consistent with those of Flambaum and Murray |12| , who 
extract from the anapole moment an f n about twice the DDH best value, f^ DHb - v - ^ 4.6, 
and point out that theory can accommodate this. (The DDH reasonable range is 0-11.4, in 
units of 1CT 7 .) However, this ignores P 7 ( 18 F), a measurement that has been performed by 
five groups. The resulting constraint is almost devoid of theoretical uncertainty 



- 0.6 < A- 0.11/4 -0.19/ii< 1.2. 



(10) 



Allowing h l p and to vary throughout their DDH reasonable ranges, one finds —1.0 < 
f n < 1.1, clearly ruling out f n ~ 9. Fig. 1 illustrates this, as well as the additional tension 
between Cs, p + a, and v4 7 ( 19 F). 

In summary there appears to be a puzzle to sort out. One hopeful development involves 
recent discussions of direct atomic measurements of anapole moments through El/ Ml in- 



terference in hyperfine transitions [19], rather than by indirect differentials. Such progress 



S 



is important because our understanding of V(e)-A(N) interactions also affects the interpre- 
tation of experiments like SAMPLE [20|, where a similar discrepancy between theory and 
experiment exists. 

This work was supported in part by the US Department of Energy and by the National 
Science Foundation. 
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